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ABSTRACT 

We investigate the angular two-point correlation function of temperature in the 
WMAP maps. Updating and extending earlier results, we confirm the lack of cor- 
relations outside the Galaxy on angular scales greater than about 60 degrees at a 
level that would occur in 0.025 per cent of realizations of the concordance model. This 
represents a dramatic increase in significance from the original observations by the 
COBE-DMR and a marked increase in significance from the first-year WMAP maps. 
Given the rest of the reported angular power spectrum C'i, the lack of large-angle 
correlations that one infers outside the plane of the Galaxy requires covariance among 
the Ci up to £ — 5. Alternately, it requires both the unusually small (5 per cent of 
realizations) full-sky large-angle correlations, and an unusual coincidence of alignment 
of the Galaxy with the pattern of cosmological fluctuations (less than 2 per cent of 
those 5 per cent). We argue that unless there is some undiscovered systematic error in 
their collection or reduction, the data point towards a violation of statistical isotropy. 
The near-vanishing of the large-angle correlations in the cut-sky maps, together with 
their disagreement with results inferred from full-sky maps, remain open problems, 
and are very difficult to understand within the concordance model. 
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1 INTRODUCTION 

Over a decade ago, the Cosmic Background Explorer Dif- 
ferential Microwave Radiometer (COBE-DMR) first re- 
ported a lack of large-angle correlations in the two-point 
angular-correlation f unction, C(8), of the c osmic microwave 
background (CMB) jHinshaw et al.1 Il99fj ). This was con- 
firmed by the Wilkinson Microwave Anisotropy Probe 
(WM AP) team in thei r analysis of their first year of 
data JSpergel et alj|2003l) and by us in the WMAP three- 
year data (|Copi et al.ll2007l) . Tho se findings have since bee n 
confirmed by iHaiianl (|2007l ) and iBunn fe Bourdon! (|2008l ). 
Here, we present a more detailed analysis of the three-year 
and (for the first time) of the five-year WMAP data, con- 
firming and strengthening our previous results. 

There is a common misconception that this lack of an- 
gular correlations is equivalent to the low quadrupole in the 
two-point angular power spectrum, which on its own does 
not have sufficient statistical significance to challenge the 
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canonical paradigm. It is typically assumed both that the 
angular power spectrum, Ci , and the two point angular cor- 
relation function, C(9), contain the same information; and 
that consequently studying one is as good as studying both. 

Actually, the exact informational equivalence between 
Ce and C(6) holds only when the full sky is observed. Statis- 
tically they are equivalent only when the sky is statistically 
isotropic. But even if Ci and C(8) did contain the same infor- 
mation, we are well aware that transforming between differ- 
ent representations of the same information — a time series 
and its Fourier transform for example — may make a real 
signal in the data easier or harder to detect. The Doppler 
peaks of the CMB, so clearly visible in the Ci representation 
are quite invisible in the two-point correlation function. 

In this work we demonstrate that outside the region 
of the sky dominated by our Galaxy, both of the CMB- 
dominated microwave bands — V and W — as well as the 
Independent Linear Combinations (ILC) map synthesized 
from them as the best map of the CMB, possess above 60° 
a level of two point angular correlation so low that no more 
than 0.025 per cent of random realizations of the best-fit 
ACDM model would have less. Indeed, above 60°, C(6) is 



© 2008 RAS 



2 C.J. Copi, D. Huterer, D.J. Schwarz and CD. Starkman 



almost entirely due to correlations involving points inside 
the Galaxy. 

This level of statistical unlikelihood (0(1O -4 )) should 
be contrasted with what could be inferred from COBE 
(O(10~ 2 )). This is a strong argument against the criticism 
that its identification as an anomaly is a posteriori. It may 
have been a posteriori for COBE, but its reidentification 
in WMAP at dramatically increased statistical significance 
is precisely how one goes about confirming that anomalies 
are actually present rather than statistical accidents of an 
observation. 

While the full-sky C(9) itself has unexpectedly low 
large-angle correlations (occurring only in 5 per cent of ran- 
dom realizations of the concordance model), we find that 
what little correlation it does have is effectively "hidden" 
behind the Galaxy. In fact, we find that a random rotation 
of the Galactic cut is as successful in masking the power 
only 2 per cent of the time, in agreement with the previous 
claim that the little correlation above 60° stems solely from 
two specific regions within the G alactic cut covering just 9 
per cent of the sky (IHaiianll2007l ). This further underlies the 
striking lack of power outside the Galactic cut, and calls into 
question cosmological uses of full-sky maps. 

Finally, we demonstrate that the absence of large-angle 
correlations is emphatically not a matter just of the low 
quadrupole. Rather, given the other measured multipoles, 
obtaining this little large-angle correlation on the cut sky 
(i.e. the part outside the Galaxy) requires carefully tuning 
C2, C3, C4, and C5. There is also a strong indication that 
it is not enough to find a model in which the theoretical 
Ce yield a very small correlation function on large angular 
scales. This is because, even if the theoretical Ce were to be 
set equal to those that are inferred from the cut-sky C(9) — 
so that the expected C(9) nearly vanished above 60° — an 
actual realization of Gaussian-random statistically indepen- 
dent ai m with these Ce) would yield different observed Ce 
because of cosmic variance. C(9 > 60°) would then not be 
nearly so close to zero. Thus getting C(9 > 60°) to vanish as 
it does seems to require covariance among the low-^ Ce, and 
thus among ae m of different i. This is in contradiction to the 
predictions of standard inflationary cosmological theory. 

One is therefore placed between a rock and a hard 
place. If the WMAP ILC is a reliable reconstruction 
of the full-sky CMB, then t here i s overwhelming ev i- 
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ments between (at least) the quadrupole and octopole and 
between these multipoles and the geometry of the Solar Sys- 
tem — a violation of statistical isotropy that happens by 
random chance in far less than 0.025 per cent of random re- 
alizations of the standard cosmology. If, on the other hand, 
the part of the ILC (and band maps) inside the Galaxy 
are unreliable as measurements of the true CMB, then the 
alignment of low-£ multipoles cannot be readily tested, but 
the magnitude of the two-point angular correlation function 
on large angular scales outside the Galaxy is smaller than 
would be seen in all but a few of every 10,000 realizations. 

We can only conclude that (i) we don't live in a standard 
ACDM Universe with a standard inflationary early history; 
(ii) we live in an extremely anomalous realization of that 



cosmology; (iii) there is a major error in the observations 
of both COBE and WMAP; or (iv) there is a major error 
in the reduction to maps performed by both COBE and 
WMAP. Whichever of these is correct, inferences from the 
large-angle data about precise values of the parameters of 
the standard cosmological model should be regarded with 
particular skepticism. 

In the remainder of this paper we provide the detailed 
evidence and reasoning behind the above statements. 



2 ANGULAR CORRELATION FUNCTION: 
PRELIMINARIES 

The two-point correlation function of the observed CMB 
temperature fluctuations 



C{9) = Tfa)T{fr) 9 , 



(1) 



where the over-bar represents an average over all pairs of 
points on the sky (or at least that portion of the sky be- 
ing analyzed) that are separated by an angle 9. On the one 
hand, we are interested in this quantity as a partial charac- 
terization of the observations. On the other hand, we regard 
it as an (unbiased) estimator of the ensemble average of the 
same quantity — where the ensemble is of realizations of the 
sky in a particular model cosmology. 

It is commonly thought that C(9) contains the same 
information as the angular power spectrum, 
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(Here ag m are the coefficients of a spherical harmonic decom- 
position of the temperature fluctuations on the sky.) This is 
because, for a full sky, 

00 

C{9) = ^^(2£+l)C^(cos6»). (3) 

e=o 

Again the Ct are regarded as most interesting to us as un- 
biased estimators of the ensemble average of |a^ m | 2 . Fur- 
thermore, the standard inflationary model predicts that the 
Universe is statistically isotropic, so that the ensemble av- 
erage of pairs of aim are independent: 
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Theoretically, the Ce therefore encode all of the information 
from the sky that has cosmological significance. 

Actually, Ce and C(9) only contain precisely the same 
information for full-sky data. Their analogues in the ensem- 
ble are informationally equivalent only if, in the ensemble, 
the sky is statistically isotropic. This suggests that by mea- 
suring both C(0) and Ce we can probe the correctness of the 
assumption of statistical isotropy of the Universe. Statistical 
isotropy is a fundamental prediction of generic inflationary 
models. 

More importantly, it is well known that, while a func- 
tion and its Fourier transform posses exactly the same infor- 
mation differently organized, in different circumstances one 
or the other may more clearly show an interesting feature. 
For example, a sharp delta-function spike in a time series 
will merely contribute equally to all modes of the associated 
Fourier series. It is precisely the same with the two-point 
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Figure 1. The two-point angular correlation function from the 
WMAP 5 year results. Plotted are C(8) for maps with Doppler 
quadrupole subtracted. The V (dashed-dotted-dotted line), W 
(dashed-dashed-dotted line), ILC (KQ75, dashed line) have had 
the KQ75 mask applied. The full-sky ILC result (solid line) is 
also shown. Also plotted are C(8) from the WMAP maximum 
likelihood Ci (dotted-dashed line), the WMAP pseudo-CV (dotted 
line) and the best-fit ACDM Of. The shaded region is the one 
sigma cosmic variance bound on the standard ACDM theory. 



angular-correlation function and its Legendre-transform, the 
angular power spectrum. Thus, while our theory may sug- 
gest to us that it is easier to analyze the angular power spec- 
trum, prudence demands that we also consider the proper- 
ties of the angular correlation function, all the more so since 
our actual measurements are done in "angle-space" not in 
"^-space" . 

In order to highlight these differences, we use the cal- 
ligraphic symbol, C, for objects operationally defined in 
"angle-space" and the symbol, C, for quantities in "^-space"; 
e.g. the Legendre transform of the two-point correlation 
function |T]) is 



2tt 



P e (cos8)C(9)d(cos8). 
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Note that Ce can be negative — in contrast to the angular 
power spectrum Ce as defined in @. 

The angular correlat ion functions in this work have been 
calculated using SpICE (|Szapudi et alj|200ll ) at NSIDE=512 
for data maps and at NSIDE=64 for the Monte Carlo studies. 
The map average has been subtracted in all cases. The re- 
sults are shown in Fig. [T]for four different maps — the ILC 
map, which covers the full sky, and KQ75 cut-sky versions 
of the ILC, the V-band, and the W-band. In the same fig- 
ure, we have plotted the Legendre transform of the angular 
power spectrum (cf. equation [3J calculated using both the 
pseudo-CV method (essentially equation [2} applied by the 
WMAP team in their first-year analysis, and the maximum 
likelihood estimates of the angular power spectrum as used 
by WMAP in the third- and fifth-year analysis. Finally we 
have plotted the expected C(9) for the best-fit ACDM, and, 
in blue, the one-sigma cosmic- variance band around the best 
fit. 

Three striking observations should be made about C(0): 



(i) None of the observational angular correlation func- 
tions visually match the expectations from the theoretical 
model. 

(ii) All of the cut-sky map curves are very similar to each 
other, and they are also very similar to the Legendre trans- 
form of the pseudo-CV estimate of the angular power spec- 
trum. Meanwhile the full-sky ILC C(0) and the Legendre 
transform of the MLE of the Ce agree well with each other, 
but not with any of the others. 

(iii) The most striking feature of the cut-sky (and pseudo- 
Ce) C(9), is that all of them are very nearly zero above about 
60°, except for some anti-correlation near 180°. This is also 
true for the full-sky curves, but less so. 

In order to be more quantitative about these observa- 
tions, we must adopt some statistic that measures large- 
angle correlations. This means that we must identify some 
norm that measures the difference between two functions 
over a range of angles. Different choices of norm, or differ- 
ent choices for the angular range, will give slightly different 
numerical results for the improbability of the above obser- 
vations; however, as we shall see, the observations are so 
unlikely that we can be confident that reasonable choices of 
the norm lead to similar results. 

In their analysis of th e first year data, the WMAP team 
defined the S1/2 statistic (|Spergel et alJuOOj ) 
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While the choice of 1/2 as the upper limit of the integral, and 
the particular choice of a square norm were a posteriori, they 
are neither optimized nor particularly special. Moreover, the 
absence of large-angle correlations was noted by the COBE 
team (though without definition of a particular statistical 
measure), and the choice of ~ 60° is clearly s uggested by 
the COBE-DMR4 results (jHinshaw et alJll996h . 

The calculation of S%/2 by direct use of is susceptible 
to noise in C(9). To avoid this we calculate Si/ 2 directly from 
Ct as 



S1/2 = CiX^e (1/2)0(1 



(7) 



The calculation of Te,e'(x) is described in Appendix lAl The 
Ce smooth over C(6) as defined in Eq. {5). 

We can use S1/2 to characterize the likelihood of the 
observed correlation function. For the COBE-DMR data 
l|Hinshaw et all 1996), there are relatively large error bars on 
C(9), which are consistent with a wide range of S1/2 ranging 
from under 1000 ()J-K) 4 to approximately 6000 (fiK) 4 . But to 
understand the significance of these values, we must compare 
them to those obtained from random realizations of the sky 
in the concordance ACDM model with the best-fit parame- 
ters. For this comparison, we generated maps based on the 
WMAP five-year ACDM MCMC parameter chain. There are 
20,401 sets of parameters in this chai n. We computed t he Ce 
for these parameter sets using CAMB (|Lewis et all 120001 ). For 
the Ce corresponding to each set of parameters, we gener- 
ated a number of random maps (i.e. maps with ae m drawn 
from Gaussian distributions with zero mean and variance 
Ce) based on the weight assigned to each WMAP MCMC 
parameter set. This produced a total of 99,997 maps at 
NSIDE=64. From the distribution of Si/ 2 values generated 
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Table 1. The Ce calculated from C(8) for the various data maps. The WMAP (pseudo and reported 
MLE) and best-fit theory C< are included for reference in the bottom five rows. 
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we calculated the probability (p-value) of randomly attain- 
ing a S1/2 as low as those we found. For COBE-DMR the 
maximum value of S1/2 of 6000 (/1K) 4 corresponds to a 3 per 
cent chance of obtaining this little angular correlation in a 
random realization of the concordance model. 



3 BASIC RESULTS 

Table [1] lists (columns 2 and 3) the value of S1/2 and its p- 
value among the sample of 99,997 WMAP MCMC maps for 
the three-year and five-year maps related to those plotted 
in Fig. [1] These include the three and five-year V band (V3 
and V5) and W band (W3 and W5) cut-sky maps. Also the 
ILC map in three-year and five-year versions, both full and 
cut sky. We use the kpO cut for three-year maps and the 
KQ75 cut in the case of the five-year maps. The five-year 
cut-sky maps are presented both as measured and corrected 
for the contribution of t he Doppler quadrupole (DQ), see 
e.g. ISchwarz et all (|2004h . (All maps are corrected by the 
WMAP team for the Doppler dipole.) 

The Legendre transform of C(9) gives us Ce, and these 
values are also listed in Table [1] (columns 4-7) for I = 2-5. 
Also included in the table, in the bottom five rows, are the 
£ = 2-5 values of the WMAP three-year pseudo-C^, the 
WMAP three- year Ce as extracted by the WMAP team us- 
ing a maximum likelihood analysis (the reported values of 
the Ce), the reported five-year values of the Ce, and the the- 
oretical Ce computed using the best-fit parameter values as 
reported both in the three-year and the five-year WMAP 
analysis. Finally and importantly, the table also shows the 
values of 5*1/2 an d their p- values computed from the Legen- 
dre transform of these angular power spectra. 

The three cut-sky maps, V, W, and ILC, whether three- 
year or five-year, are all in good agreement with each other. 
They all have very low values of S1/2 — almost two orders- 
of-magnitude below the predictions of the theory. In both 



the three and five-year ILC outside the Galaxy, the prob- 
ability that such low values could happen by chance is ex- 
traordinarily low — only 0.025 per cent. This low probabil- 
ity is e ntirely consistent wit h the original finding of COBE- 
DMR i|Hinshaw et al.1 Il996l ) described above, however the 
error bars on C(0) (and hence on S1/2) have declined sub- 
stantially, with the WMAP value of S1/2 being at the abso- 
lute lowest end of what was consistent with COBE, and with 
much smaller error bars. This dramatic decline in the error 
bars, while honing in on the very low end of the COBE-DMR 
range, is exactly what one would expect from the absence 
of large-angle correlations in the CMB sky, and not at all 
what one would expect if the low S1/2 in COBE-DMR (and 
in WMAP) was merely a statistical fluctuation in the mea- 
surement. 

We also consider the case where there is exactly zero 
large scale angular correlations. That is, we set C(8 ^ 60°) = 
and extract the Ce as a Legendre transfornQ for the ILC 
(kpO) map. This "theory" produces low-^! Ce of approxi- 
mately the same value as for the C(8) from the cut-sky maps. 
On the one hand this is consistent with the statement that 
there is little correlation on large angular scales and thus 
the Ce for \ow-£ are dominated by small angular scales. On 
the other hand, this shows that the data is consistent with 
there being no correlations on large angular scales. 

We note that it is difficult to enforce C(9) =0 in the con- 
text of a statistically isotropic model. Even if a model were 
found that predicted the observed Ce as the expected means 
of the \ae m \ 2 (as in equation fj|, any actual realization of the 



1 We note that setting C(6 60°) = introduces a small 
monopole into the power spectrum. This can be corrected by 
subtracting it out, changing the 8 C = 60° to a value such that 
J C{6) smOdO = 0, etc. Without an underlying theory to de- 
scribe this case it isn't clear how to best correct for this monopole. 
Regardless, the Ce we extract are not very sensitive to the method 
we use for removing the monopole. 
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Universe would produce Ci that were substantially different. 
Indeed, we have found that > 97 per cent of realizations of 
such a Universe would have values of S1/2 greater than the 
observed value (see section [5~Tj) . 

The results from the full-sky ILC map, also show low 
values of S1/2; however, with less remarkable p- values of 5 
per cent in contrast to 0.02-0.04 per cent for the various 
cut-sky maps. Similarly, the full sky ILC maps have larger 
quadrupoles than the cut-sky maps (though still lower than 
expected from theory), and octopoles consistent with theory. 
These full-sky maps are in good agreement with the WMAP- 
reported MLE Ce- Meanwhile the pseudo-C^ based on the 
kp2 sky-cut (which cuts less of the sky than the kpO cut) 
are intermediate between the kpO cut-sky map results and 
the full-sky results. 

Thus the full-sky results seem inconsistent with cut- 
sky results and they appear inconsistent in a manner that 
implies that most of the large-angle correlations in recon- 
structed sky maps are inside the part of the sky that is con- 
taminated by the Galaxy. 



4 MISSING POWER OR UNFORTUNATE 
ALIGNMENT? 

An important question to consider is whether the extremely 
low large-angle correlations in the cut-sky WMAP maps are 
a general result of cutting the maps or is specific to the 
orientation of the cut. That is, should we expect a loss of 
large-angle correlations in a cut-sky map or is the alignment 
of the cut with the Galaxy important. To address this ques- 
tion the full-sky five-year ILC map was randomly rotated 
(that is, set its north pole in a random direction and draw 
its azimuthal angle from a uniform distribution) 300,000 
times. For each random rotation we masked the map with 
the Galactic KQ75 mask, which is now, therefore, randomly 
oriented relative to the original map and re-computed the 
quadrupole, octopole and 51/2 statistic. 

The analysis shows that the true cut-sky quadrupole 
and octopole are not terribly unusual compared to those in- 
ferred from the rotated-then-cut (RTC) maps. In 7.6 per 
cent of these RTC maps the quadrupole is smaller than 
that of the ILC with the originally placed cut, while 2.5 
per cent have a smaller octopole. Therefore, if we looked at 
the quadrupole and octopole alone we would conclude that 
an arbitrarily oriented mask is only moderately unlikely to 
produce low-^ power in the cut-sky ILC. Conversely, the par- 
ticular alignment of the Galaxy with the part of the sky on 
which the low-£ power is concentrated is only moderately 
important. 

On the other hand, in the RTC maps 5*1/2 = 
11900 (mK) 4 with variance 7300 (fiK) , a very high value 
relative to the original cut ILC map (1152 (/iK) 4 , see Ta- 
ble QJ. Only 2 per cent of these rotated maps have S1/2 
lower than the ILC with the original cut. (Recall that 
5i/2 — 8583 (/iK) 4 for the full-sky ILC is already low, with 
a p-value of only about 5 per cent.) Thus it is quite unlikely 
for an arbitrary cut to suppress the large-angle correlations 
to the extent observed in the cut-sky ILC map. Conversely, 
it is quite likely that the observed absence of large-angle cor- 
relations outside the KQ75 cut is due to the alignment of the 
Galaxy with the regions on the sky where such correlations 
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Figure 2. The two-point angular correlation function from the 
WMAP 5 year results. Plotted are C(9) for the ILC calculated 
separately the part of the sky outside the KQ75 cut (dashed line), 
inside the KQ75 cut (dotted line), and on the part of the sky with 
at least on point inside the KQ75 cut (dotted-dashed line). For 
better comparison to the full-sky C(9) (solid line), the partial-sky 
C(0) have been scaled by the fraction of the sky over which they 
are calculated. 

are maximi zed. Th i s resu lt is in good agreement with the 
result from iHaiianl (|2007l ) that the little correlation above 
60° stems from two specific regions within the Galactic cut 
covering just 9 per cent of the sky. 

It appears that our microwave background sky has 
anomalously low angular correlations everywhere outside 
the Galactic mask, but not within. In Fig. [2] we plot C{6) 
for the WMAP 5 year ILC map calculated separately on the 
part of the sky outside the KQ75 cut, inside the KQ75 cut, 
and on the part of the sky with at least one point inside the 
KQ75 cut. For better comparison to the full-sky C(6) (also 
plotted), the partial-sky C{9) have been scaled by the frac- 
tion of the sky over which they are calculated. This shows 
that the full-sky C(9) is very close to C(8) calculated from 
the masked region. Meanwhile C(8) calculated outside the 
Galactic mask is similar to neither, and much closer to zero 
in magnitude. 

Figure shows two other interesting peculiarities of the 
measured angular correlation function. First, the full-sky 
C(9) is particularly closely mimicked by the C(9) computed 
so that at least one of the points is inside the masked re- 
gion; the agreement between the two at large angles (above 
approximately 60°) is nearly perfect. Moreover, all four cor- 
relation functions shown in the figure vanish at nearly the 
same angle, 6 ~ 90°. While at this time we do not under- 
stand the origin or significance of these two effects, we wish 
to point them out because it is possible that they will be 
useful for successful theoretical or systematic explanations 
of the vanishing correlation function. 

The evidence we present strongly suggests that the full- 
sky ILC map does not represent a statistically isotropic mi- 
crowave sky. If the region outside the cut is a reliable repre- 
sentation of the CMB then we should focus on the angular 
correlation for cut skies. As shown above this leads to a p- 
value of 0.025 per cent for the standard ACDM model (see 
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Table [TJ. Furthermore, the WMAP reported MLE C e are in 
good agreement with the full-sky Ce and C(9), but not with 
their cut-sky equivalents, whereas, the cut-sky Ce and C(9) 
are in good agreement with the pseudo-Ce (see Table\l\ and 
figure QP. This casts doubt on the validity of the reported 
low-i Ce- 



5 SO, IS THE LARGE SCALE CMB 
ANOMALOUS OR NOT? 

It has been suggested that there is not hing particularly 
anomalous abou t the l arge-angle CMB (lEfstathioul |2004| ; 
iGaztanaga et al.l |20o3 : ISlosar etail 120041 ). The argument 
goes something like this: (a) the two-point angular correla- 
tion function C(9) and the angular power spectrum Ce con- 
tain the same information; (b) not only does theory tell us 
that that the Ce are the "relevant" variables, but, since they 
are discrete and finite in number, we can apply standard sta- 
tistical techniques to compare observations with theoretical 
predictions; when we do so, we find that only C2 is far be- 
low the expected value, but still at a level that happens by 
chance 10 per cent of the time in the concordance model. 

We have already pointed out above that (a) contains 
two fallacies. First, C{9) and Ce are equivalent only on a 
full sky. Second, formal equivalence is not the only question, 
signals are often more visible in one representation of the 
data than in a different, though formally equivalent one. If 
this were not so, then there would be no need for Fourier 
analysis. Nor would we need to perform great music — we 
could simply read the score. 

Point (b) is correct if the Inflationary Cold Dark Mat- 
ter (ACDM) cosmological model is true. This model tells us 
that the spherical harmonic coefficients ae m are independent 
Gaussian random variables, and that the sky is statistically 
isotropic, so that the off-diagonal covariances of the estima- 
tor Ce (as defined in Eq. vanishes in linear theory. This 
vastly simplifies statistical analysis of the CMB in the con- 
text of ACDM. However, in advancing the case for a particu- 
lar cosmological model we are required not just to determine 
the best-fit parameters of the model but to test the assump- 
tions and other predictions of the model. This includes the 
prediction of statistical isotropy, and consequently that the 
aim are independent of one another. 

There is already considerable evidence that if one an- 
alyzes the full-sky ILC that one finds difficult-to-explain 
deviations from statistical isotropy, such as the align- 
ment of the octopole and quadrupole with each other 
and with the geometry of the S ol ar System (fo r exam - 
ple Ide Oliveira-Costa et al.1 ll2004h; lEriksen et al.1 l|2004h ; 
ISchwarz et all (|2004h ; lLand fc Magueiiol (|2005aF n. These 
analyses require full-sky data for any statistical power, and 
so, in particular, might be explained by Galactic foreground 
contamination (though it would be an odd coincidence for 
Galactic contamination to cause alignment with the Solar 
System). The calculation of C(8), as we have seen, can be 
done on a cut sky. 



5.1 Are just the low-£ C'e incorrect? 

Perhaps the standard assumptions of statistical isotropy and 
Gaussianity are correct and only the low-^ Ce are incorrectly 



predicted by the standard model. If this were the case then 
the ae m could still be Gaussian random variables and some 
new physics would be needed to explain the low-^ Ce, e.g., by 
giving up the scale invariance of the primordial power spec- 
trum, which could be caused by a feature in the inflationary 
potential. 

To study this possibility we replaced C2 through C20 in 
the best fit ACDM model with the values extracted from the 
cut-sky ILC five-year map. From these Ce 200,000 random 
maps were created, masked, and S1/2 computed. Under the 
assumptions of Gaussianity and statistical isotropy of these 
Ce only 3 per cent of the generated maps had S1/2 less than 
1152 (^K) 4 (the cut-sky ILC5 value from Table [lj. Thus 
even if the Ce are set to the specific values that produce 
such a low S1/2 a Gaussian random, statistically isotropic 
realization is 97 per cent unlikely to produce the observed 
sky. Again this shows that either (i) the \ow-£ Ce are corre- 
lated, thus breaking statistical isotropy, or (ii) our Universe 
is a 97 per cent unlikely realization of an alternative model 
that deviates from the standard one as required to produce 
the Iow-I Ce- 

5.2 Statistics of C{6) 

Once we have decided to calculate C(8), we are forced to ask 
how best to analyze it statistically. One option would be to 
compare the C(9) inferred from a particular map to the C(9) 
one expects from theory. Thus one would define 



C th (8) = V(2*+l)Cf P*(cob0), 

47T t—^l 



(8) 



for a particular set of parameters (say the best-fit values) of 
the concordance model. This is what is plotted in figure [T] 
as the C{6) of the best-fit ACDM model. 

One would next define some functional norm and com- 
pute 

ohs - ACDM = \\c(e)-c th (e)\\ (9) 



N 



where we imply a suitable average over a range of 9 on the 
right-hand side. This norm could serve as a statistic to com- 
pare the two-point correlation function inferred from the 
data, or some subset of the data to the theory. The shaded 
band around C th (9) in figure Q] (cosmic variance) reflects 
this notion that one somehow expects the inferred C{9) to 
lie inside this band. 

The statistic originally suggested by the WMAP team 
for comparing observations of large-angle correlations to the- 
ory, S1/2, does not fall into the above class of statistics. This 
is because it captures that what is strange about the inferred 
angular correlation function C(9) is not that it is different 
than theory for 9 > 60°, but rather that it is so close to 
zero. Thus S1/2 is designed to test an alternative simple hy- 
pothesis — that there are no correlations above 60°. In the 
language of equation © S1/2 is in the class of statistics 



N obs - ZCI ° = \ \C(9) — 0||. 



(10) 



There is another lesson to be learned from the preced- 
ing results. Cosmological inflation predicts that there are 
fluctuations on all scales, whereas many alternative mod- 
els of structure formation, like cosmic defects, would pre- 
dict the absence of fluctuations on super-horizon scales. By 
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looking at scales above 1 degree on the sky the inflation- 
ary prediction is tested at the time of photon decoupling, 
and by looking at the largest angular scales, we can test it 
in the more recent Universe since the physical Hubble scale 
Th{z) = 1/H(z) is observed at the angle 8 = r H (z)/d a (z) 
and angular distance d a (z) = [1/(1 + 2)] L [1/H(z')]dz' . For 
the best-fit parameters of the concordance model, the lack 
of correlations at larger 60 degrees means, that scales that 
crossed into the Hubble radius below a redshift ~ 1.5 are 

uncorrelated. 

Instead of S1/2, iHaiianl (|2007T ) advocates the use of a 
covariance- weighted integral over C(8), 

d(cos0) / d(cos 8' )C (0) F' 1 (6, 6' )C(0'), (11) 

where 

F(8, 8') ee {{C{6) - <C(0)»(C(0') - <C(0')») , (12) 

and (• • •) represents an ensemble average, i.e. an average 
over realizations of the underlying ACDM model. As Ha- 
jian notes, .4(1/2) = S1/2 in the limit of uncorrelated C(8). 
However, just because C(8) and C(8') axe correlated in the 
standard theory does not make 4(1/2) a more correct statis- 
tic than Si/2- For tests against the standard theory the A(x) 
statistic provides another statistic; one that accounts for the 
theory correlations. However, as discussed above, it has re- 
peatedly been shown that there are correlations among the 
low-£ multipole moments (and multipole vectors) of the full 
sky that are not consistent with the standard theory. In this 
case it is not possible to compute F(9,0') because the en- 
semble over which to average is unknown. T herefore , while 
it is somewhat reassuring t hat by using A(i ) IHaiianl (|2007l ) 
confirms our earlier result (|Copi et al.ll2007l ) showing C{8) 
for cut skies violates the fundamental model assumption of 
statistical isotropy, it is not clear that any strong inference 
should be drawn from differences in statistical significance 
between results for 4(1/2) and Si/2- 

Moreover, we again emphasize that what is anoma- 
lous about the observed large-angle correlations is not how 
poorly they match the theory, but rather how well they 
agree with the very simple alternative phenomenological 
hypothesis that there are no large-angle correlations - 
C(8 > 60°) = 0. The construction of jV obs - ACDM , m i g ht 
indeed benefit from an attem pt to remove e xpected corre- 
lations through F(8,0'), as in IHaiianl i|2007h ; however, the 
theoretical model against which jv obs_zcro compares the ob- 
servations has F(8, 8') = for the relevant 8. 



5.3 Minimizing S1/2 

Once we have understood that what is anomalous about 
C(8) is how close it is to zero, we can understand that what 
is strange about the low-^ Ce is not just how low C2 is, but 
also how the various Ce are correlated with each other. 

We now probe the sensitivity of S1/2 to ranges of I. 
Given that small angles can affect low-^ results, it is also the 
case that higher I can affect the larger angles. One way to see 
this is to determine how the Ct for low-^ must be adjusted 
to attain a low SV/2. This is not done by setting some range 
of Ct to zero. Instead, given a set of Ce for I > f max ,timc we 
can find the values of Ct for 2 ^ I ^ ^max.tuno that minimize 



S1/2 by regarding S1/2 as a function of Ct using equations 
© and ©. 

We consider two sets of Ct, the first from the ACDM 
theory, the second as reported by WMAP. Table[2]shows the 
minimum S1/2 we find for each value of 4nax,tune- In the ta- 
ble, we provide the values for the best-fit ACDM model and 
the reported WMAP Cg. We also provide the 95 per cent 
confidence ranges based on the WMAP MCMC parameter 
set chain where the minimum S1/2 was calculated indepen- 
dently for each model. 

To attain S1/2 ^ 1152 (fJ-K) 4 (the value found in the 
masked ILC map, see Table [TJ from the best-fit theory re- 
quires tuning both C 2 and C3. Thus even the theory requires 
more fine-tuning than just the quadrupole to be low in or- 
der to be consistent with observations. The minimum in Ta- 
bled] was attained for 6C* 2 /27r = 149 (^K) 2 and 12C 3 /2tt = 
473 (/1K) 2 . In general for the theory we need to tune at least 
up to f m „,tunc = 3 and can almost always find a low S^/2 if 
we tune up to W,tunc = 4. 

For the WMAP Ct even more tuning is required. Note 
that the WMAP C2 is already approximately tuned to 
produce the minimum S\/2 given the rest of the Ct for 
t > 2 (that is, from Table [1] we note that 8583 (^K) 4 w 
8290 (fiK) 4 ). To attain the low Si/ 2 to match the cut-sky 
ILC requires tuning of values of Ct up to ^max.tuno = 5. 

Table [2] further shows that the minimum 5*1/2, that can 
be achieved by optimizing the \ow-£ Ct fall off much more 
slowly in the WMAP Ce than in the theory. By ^ m ax,tunc = 
8 the minimum WMAP S1/2 is two orders of magnitude 
larger than can be attained from the theory. This strongly 
suggests that important correlations exist in the data for 
I 8 that do not exist in the theory. These correlations 
cannot be canceled by tuning the lower £ behavior. 

Therefore, we conclude that a given behavior of C(8) 
on large scales is not in unique relation to a behavior of 
Ce at low-^. The former quantity receives significant con- 
tributions from Ce at high I as well; the converse is also 
true. Given the extremely puzzling near-vanishing power in 
C(8 > 60°), and given that the qua drupole and octopole a re 
not unusually low (as shown in e.g. lO'Dwver et al.l (| 20041) 1. 
we argue that any theoretical or observational explanation 
of the "low power at large scales" should concentrate on the 
quantity C(8). 



6 CONCLUSIONS 

In this paper we have studied the angular correlation func- 
tion in WMAP three- and five-year maps. We have clarified 
the relation between various definitions of the angular corre- 
lation func t ion, a nd revisited our previous calculation from 
ICopi et all i|2007l ') in more detail. We confirmed that power 
on large angular scales — greater than about 60 degrees — is 
anomalously low, at 99.975 per cent CL (see Table [1]). The 
measured angular correlation function thus disagrees with 
the ACDM theory, but, more significantly, it is consistent 
with a simple phenomenological "theory" — C(8 60°) = 0. 
The significance of this disagreement (as measured by the 
probability of the value of S1/2) has now increased by a fac- 
tor of over 100 since it was first observed in the COBE-DMR 
four- year analysis. 

We have shown that the cut-sky and full-sky large-scale 
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Table 2. S]/jj (in (/jK) 4 ) obtained by minimizing with respect to CV We show the 
statistic for the best-fit theory and WMAP, as a function of the cutoff multipole £ max ,tune 
(the minimization has been performed by varying all t in the range 2 $J I ^ £ m ax,tune and 
fixing I > ^max,tunc)- Also shown is the 95 per cent confidence region of the minimized 
S\ i2 derived from chain 1 of the WMAP MCMC parameter fit. In the bottom row, we 
remind the reader that the measured value of S1/2 outside the cut is 1152 (/iK) 4 (see 
Table [T] for more details). 



Source 




Maximum tuned multipole, £ ma x 


,tunc 
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WMAP 


8290 


2530 


2280 


800 


350 


150 130 


ILC5 (KQ75) 






1152 









angular correlations differ (see Table[T]and Figure[T]), though 
the source of these discrepancies remains unknown. This 
shows that either the Universe is not statistically isotropic 
on large angular scales or that correlations are introduced in 
reconstructing the full sky from the observations. We have 
shown that even given the unusually small full-sky angular 
correlations (95 per cent unlikely) an unusual alignment of 
the Galaxy with the CMB (2 per cent of realizations) is re- 
quired to explain the lack of correlations outside the Galac- 
tic region. We have further shown that simply adjusting the 
theoretical values of the Ci does not solve the problem if 
the sky is representative of a Gaussian random statistically 
isotropic process - the cosmic variance in the Cg is such that 
less than 3 per cent of all realizations would preserve a low 
value of Si/ 2 . 

From these results we argue that C(8) is an important 
quantity to study along with the usual angular power spec- 
trum, Ci. The typical "rule-of-thumb" that \ow-i describes 
large angular scales is not accurate. Any theoretical expla- 
nations for the "missing large-scale power" should concen- 
trate on explaining the low C(9), rather than the smallness 
of the q uadrupole and octopo le, which are not nearly as sig- 
nificant (lO'Dwver et aL| 2004). As has been pointed o ut by 
iRakic fc Schwarzl (|2007h and iBunn fe Bourdon! (|2008l V any 
possible explanation of the multipole alignments that relies 
on an additive, statistically independent contribution to the 
microwave sky on top of the primordial one, increases the 
significance of the lack of angular correlation. 

The CMB, as measured by WMAP in particular, pro- 
vides much support for our current model of the Universe. 
It also points the way toward new puzzles that may affect 
fundamental physics. On the largest angular scales the mi- 
crowave sky is inconsistent with theoretical expectations. 
These discrepancies between observations and theory re- 
mains an open problem. In the future, combining the current 
data with new information, such as new data from WMAP, 
observations from the Plan ck exp eriment, and polarization 
information (Dvor kin et al. 1 120081 ) may be key to determin- 
ing the nature of the large-scale anomalies. 
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APPENDIX A: INTEGRATING PRODUCTS OF 
LEGENDRE POLYNOMIALS 

We wish to calculate 

P m (x')P n (x')dx'. (Al) 

For the special case of x = 1 this is just the normalization 
2 

^m.ri(l) = ~ ; — r^m.ro- (-^-2) 

In + 1 

For a general x we consider two cases. When m ^ n 
Legendre's equation 

(1 - x 2 )P^(x) ~ 2xP„(x) + n(n + l)P n (x) = (A3) 
allows us to write 

P m (x)P n (x) = A [(l- x 2 )(P^P„-P;P m )] 

/ [n(n + l) -m(m- 1)]. (A4) 
Then using the relation 

(1 - x 2 )P^(x) = nP„_i(a;) - nxP n (x) (A5) 
we find 

lm,n = {mP„(x)[P m - 1 (x) - xP m (x)] 

- nP m (x) [Pn-i(x) ~ xP n (x)]} 
J {n(n + 1) - m(m + 1)} [form/n], (A6) 

When m = n we integrate (IA1[I by parts to get off 
diagonal terms <|A6[) and use the indefinite integral 

J Pn(x)dx = — ^ [Pn+l(x) - P„-l{x)] (A7) 

to derive the recursion relation 

1 n ,n{x) = {[P„+i(as) - P n -i(a:)] [P»{x) - P»- 2 (a:)] 

- (2n - l)X n+ i,„_i(a;) + (2ri + l)I n , n -i{x) 

+ (2n-l)J„_i,„_i(x)} / {2n + l}. (A8) 

Starting from Ta,o(x) = x + 1 and T\ t \(x) = (x 3 + l)/3 we 
can calculate all the diagonal terms recursively. 

With these two relations ()A6l and lA8|l we can compute 
and tabulate all required values of I m ,n for any x. 

This paper has been typeset from a TpX/ F/IpX file prepared 
by the author. 
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